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CN ■ Abstract 

We study theories of gravity with non-minimal coupling between polarized media with pole-dipole and 
| ^ | quadrupole moments and an arbitrary function of the space-time curvature scalar R and the squares of 

the Ricci and Riemann curvature tensors. We obtain the general form of the equation of motion and show 
that an induced quadrupole moment emerges as a result of the curvature tensor dependence of the function 
coupled to the matter. We derive the explicit forms of the equations of motion in the particular cases of 
coupling to a function of the curvature scalar alone, coupling to an arbitrary function of the square of the 
Riemann curvature tensor, and coupling to an arbitrary function of the Gauss-Bonnet invariant. We show 
that in these cases the extra force resulting from the non-minimal coupling can be expressed in terms of 
the induced moments. 
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■ 1 Introduction 

Various extensions of the general theory of relativity have been proposed in recent years to explain astrophysical 
observations like the accelerated expansion of the Universe and rotation curves of spiral galaxies. A widely 
known class of such proposals comprises variants of the so-called f(R) theories of gravity; see [TJ [H [3] for 
reviews. Even these theories have been themselves subject to different modifications in the hope of obtaining 
a more satisfactory picture, namely, by considering non- minimal couplings with matter [?J [SJ |H1 [7]. In the 
models investigated in these references some powers of the curvature scalar R have been coupled to the matter 
Lagrangian which have resulted in modified equations of motion. A more general case was considered in [5] in 
which an arbitrary function of the curvature scalar non-minimally coupled with the matter Lagrangian. Taking 
perfect fluids and scalar fields as matter fields, several implications of such a coupling were subsequently studied 
in [SJ [TU1 US HH EH H31 E3 EH H7J • The effects of such a coupling on the dynamics of test particles have been 
studied in |18j where it has been shown that an extra force emerges from the non-minimal coupling. The 
explicit form of the force depends on the form of the matter Lagrangian and it has been shown in [TOJ [5D] 
that for the case of the perfect fluid for which there are several equivalent Lagrangians, one may end up with 
vanishing or non- vanishing extra forces depending on the chosen form of the Lagrangian. In |21) , it has been 
argued that the Lagrangians which are equivalent for a non-coupled perfect fluid are not in fact equivalent 
when the perfect fluid couples non-minimally to curvature. The author of [22) has argued that the form of the 
coupling fixes the matter Lagrangian. 

The extra force affects the world-lines of test particles. In fact within the framework of models with non- 
minimal coupling, test particles move along non-geodesic trajectories in general [18] . Taking the effect of extra 
force into account, the motion of test pole-dipole particles, i.e. test bodies with certain microstructures, has 
been studied in [23 by using the method of multipoles [MJ [35]. In this work we consider the motion of pole- 
dipole and quadrupole particles in extended gravity models with more general coupling to matter. To this 
end we consider a class of matter fields with a general dependence on a tetrad field, its covariant derivative, a 
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velocity field, curvature tensors, and a generic set of matter fields, coupled to f(R) gravity action. In particular 
we show that when the Lagrangian depends on curvature tensor through non-minimal coupling to geometry, 
induced quadrupole moments emerges. We obtain the equations of motion for media with pole-dipole and 
quadrupole moments non-minimally coupled to an arbitrary function of the curvature scalar and as a special 
case we consider a Wcyssenhoff-type spinning fluid from which we extract the equations of motion of polc-dipolc 
particles. The method can then be easily generalized to particles with electric charges. We also investigate the 
motion of test bodies non-minimally coupled to an arbitrary function of the square of the Riemann curvature 
tensor or to an arbitrary function of the Gauss-Bonnet invariant and show that the extra forces can be expressed 
in terms of the induced quadrupole moments. The problem of induced quadrupole moments in media due to an 
applied gravitational field has been studied in [26, 27, 28, 29, 30 in the context of general relativity. The study 
of motion of quadrupole particles is well motivated by their role in gravitational radiation. Several examples of 
motion of quadrupole particles in curved space-times of general relativity may be found in |31[ 1321 133[ I34| 135] . 
The effect of multipole moments of planets on bending of light rays has been discussed in [315] , and it has been 
suggested in [37J [35J that such effects may be detected by using microarcsecond astronomical interferometers. 

In the subsequent sections we first give a brief review of the polarized media as formulated in [39 , 40 and then 
deploy it to obtain the equations of motion of polarized media with pole-dipole and quadrupole moments non- 
minimally coupled to an arbitrary function constructed out of the curvature scalar and the squares of the Ricci 
and Riemann curvature tensors. We subsequently consider three particular cases, the case of non-minimal f(R) 
gravity, the case of coupling with an arbitrary function of the Riemann curvature tensor squared, and the case 
of coupling to an arbitrary function of the Gauss-Bonnet invariant. In the first case, our investigation results in 
equations of motion of particles with pole-dipole and quadrupole moments in non- minimal f(R) gravity. In the 
second and third cases in which for simplicity no intrinsic moments (i.e. those originating from the structure 
of the body) are assumed, we show that the extra force of the non-minimal coupling can be expressed in terms 
of the induced moments. In the last section we discuss the results. Throughout the work we use the following 
(anti)symmetrization conventions: A"^ u >" = A"^"'" + A" Vfl "' , A" « />>■•' = ,r " /; ••' • y <■ n _ 
and a I! ! ' = A " IS " - .1 •" // " . 



2 Particles with pole-dipole and quadrupole moments in General 
Relativity 

It has been shown in [39] [40] that by starting from a Lagrangian with the general form e£, e£, R^ a p-y) 

in which it' 1 is a velocity field, is a tetrad field satisfying e a jl e > l /(] a b — g^, with r\ a b — diag(— 1, 1, 1, 1), 
= ii Q V a e^, and i? p Q( 3 7 = d a r^ u — dpT^ v + T^r^ — T^r^ is the Riemann curvature tensor, one can 
arrive at the following energy-momentum tensor 

T^y = + Vh^ v - + A^ v (1) 

describing a polarized fluid with pole-dipole and quadrupole moments. Here ii M is the four-velocity of a fluid 
element, p M its four- momentum, V is the pressure, 



= ^V Q (s Q (M u t/) + s^u a ), (3) 



hfiu = g^ + u^Uv, (2) 
1, 
2 

V = q^R^ 1 -3q va ^R^ -2VpV a q a (p „/ (4) 
where is the fluid spin density and 

= (5) 

is the quadrupole moment tensor, which is assumed to share the symmetries of the Riemann curvature tensor. 
The above energy-momentum tensor is symmetric provided the spin evolution equation 



W a (u a s^) = p^u" - P v u» - 4R^ a ^q^ a0 (6) 
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holds. The translational equation of motion 

u"V„p" +p"V„u" + WVvV + VV v (u»u") = --Rf^u's'* + q^s^^R^ 5 (7) 

can be obtained from V JT^ = or by variation of world-lines [20] • The so-called Frenkel condition 

u^aT = (8) 

ensures that spin remains space-like. Imposing the energy conditions may put further restrictions on the spin 
tensor. By neglecting the quadrupole moments, the equations of motion of a spinning fluid (see e.g. |41j and 
references therein) are recovered. By turning the spin off and setting the quadrupole moments equal to zero, 
the standard perfect fluid energy- momentum tensor is recovered from Eq. (TT]). 

As has been shown in [391 140) . the above procedure can be extended to Lagrangians of more general 
dependence of the form 

£ = />V£,e£,<^) (9) 

in which 

<f> A = {R\^, V M1 • • • V Mm i?V 7 , </, V Q </} (10) 

with ip A (A is a collective index) being a generic set of matter fields whose free field Lagrangians are also 
included in the total Lagrangian. In particular by taking (p A to be the electromagnetic field, the equations of 
motion for charged particles would be obtained. Thus in the special case where particles have only pole-dipole 
moments, the Mathisson-Papapetrou-Dixon |42[ and Dixon-Souriau equations [43] for charge-less and charged 
particles are recovered respectively. Multipole moments higher than the quadrupole moment are incorporated 
via dependence on covariant derivatives of the Riemann curvature tensor. 



3 Multipole particles in non-minimal extended gravity 

The action for f(R) gravity with non-minimal coupling to matter is given by [18] 

S = J + U + A/ 2 (i?)}£ ro J J—gd A x (11) 

where C m represents the matter Lagrangian, A is a coupling constant, and fi(R) and $2{R) are arbitrary 
functions of the curvature scalar R = g ai g@ s R a p 7 s- For A = this reduces to the usual f(R) gravity. As has 
been shown in |18[ . one can obtain the following relation from the equations of motion associated with the 
above action 

VJF^ = i X + ! ^ R) (<TXn - T^) V„R (12) 
where f 2 (R) = ^-g^ , and T^ v is the matter energy-momentum tensor 

TfJ ,u = 2 S{£ m i/=g) 

Here we consider a more general extension of the above action with the following form 

\R+\f + (l + \F)C^d i x (14) 

where / = f{R 7 P 7 Q), F = F(R,P,Q), and P = R^R^, R^ being the Ricci tensor, Q = R^pR^^ is 
the Kretschmann invariant, and £ is a matter Lagrangian (see [44] for a cosmological model with this type of 
action but with A = 0). In fact one may think of (1 + \F)C as a Lagrangian of the type described by Eq. (|9]). 
Various models including minimal or non- minimal f(R) gravity, modified Gauss-Bonnet gravity, the so-called 
f(R, G) models [45J 06] where Q = R 2 — 4P + Q is the Gauss-Bonnet invariant, and the model introduced in 
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[47] can be recovered as special cases by choosing appropriate forms for / and F. Functions with more general 
dependence on curvature invariants (e.g. those introduced in |48j ) may also be considered as toy models. 

The equations of motion can be obtained by varying the above action with respect to metric components. 
Defining 

= (R^ - + g^U)f R + 2R» a R v J P - V ce V (iy (/ P i? Q ^ ) ) + g^V a Vp(f P R ap ) 

+D(/ P i?n + 2R^R v aM f Q - 2V a V ■p(f Q K*Mf>) 

we have 

Qitu + H ?u _ _f g ^ = (i + \F)T^ - 2XW V (15) 

where □ = V^V, j|, ^, % , and ^ stand for f R , f P , f Q ,F R , F P , and F Q respectively, and F is 
obtained from H^ LV by substituting 

Sr -> /p -)■ £F P , f Q ->■ £F Q . 

The matter equation of motion may be obtained by taking the divergence of both sides of Eq. (|T5l) . which 
results in 

v ^ = i -T5 (Vv-FT"" - 2V„-HH • (16) 

1 + A-T 

Now if one chooses a Lagrangian of the form given by Eq. (|9]), the contributions from intrinsic pole-dipole and 
quadrupole moments can be tracked through T^ u terms in the above equation. However even if we choose 
a Lagrangian involving no pole-dipole or quadrupole moments, quadrupole moment contributions would still 
appear, through F Rl Fp, and Fq terms, in the dynamics. Thus the non-minimal coupling induces quadrupole 
moments in the medium. The appearance of these induced moments is an essential feature of the models with 
non-minimal coupling between geometry and matter. It is a novel phenomenon not observed in other modified 
theories of gravity. To proceed further, we consider some rather interesting special cases of this model. 



4 Pole-dipole and quadrupole particles in non-minimal f(R) gravity 

In this section we restrict ourselves to a particular case of the model described above where / = f(R), 
F = F(R), and C = C{u^, e£,e%, i? M a ^ 7 ); i.e. we consider particles with pole-dipole and quadrupole moments 
in non-minimal f(R) gravity. 

Inserting Eq. ((T|) into Eq. (IT51) we obtain 

u v V vP » + p»V v u v + h^VuV + VV v {u»u v ) + ^Rfuapu^ - q^s^R^ 5 = r (17) 

with 

F = (SSCm - P »u v - VK + B»„ -A» u ) (18) 

being the extra force resulting from the non-minimal coupling. There are some quadrupole contributions to 
this extra force coming from dependence on curvature. These induced moment contributions can be computed 
with the aid of the following relation 

fir 1 

= {1 + XF)-^— + -g^g^\F R C (19) 

which is obtained from Eq. (JS|) by replacing C with (1 + XF)C. In this equation the first term in the right-hand 
side is proportional to the medium intrinsic quadrupole moment while the second one represents the quadrupole 
moment induced as a result of the non-minimal coupling. One may use the above type of equations to obtain 
a rough estimate of the magnitude of these induced moments. In the case of the model considered in this 
section the magnitude of the induced quadrupole moment relative to the intrinsic one is of the same order as 

TTXf' 
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To obtain the equations of motion of particles with pole-dipole and quadrupole moments, in Eq. (|T7| we 

set 

V = 0, p" = nP», = nS^, q\ M = nQ\ M 

and use V a (nu a ) = 0, where n is the particle number density and P' 1 , S^ u , and Q tl a p 1 are the single-particle 
momentum, spin, and quadrupole moment respectively. This results in 

tt'VtfP" = /" - ^uafiWS^ + Q Q p jS VR a ^ S . (20) 

If we neglect the quadrupole moments, this equation is the non-minimal f(R) gravity counterpart of the 
Mathisson-Papapetrou-Dixon equation and agrees with the results of [23]. By turning the spin off, the non- 
minimal coupling of / (P)-perfect fluid considered in [18] is achieved. For pole-dipole particles with electric 
charges one can add the Maxwell energy-momentum tensor to the right-hand side of Eq. (fTJ) and proceed as 
above. This would result in the non-minimal f(R) gravity counterpart of the Dixon-Souriau equations. 

A particular case of the model described above is the case where F(R) = R. If we also assume, for 
simplicity, that V — and the medium has no intrinsic pole-dipole or quadrupole moments, Eq. ([T7| together 
with Eqs. (TH)) and (JTSJ) give 

V ^ V) = rr^G^ Q ^~^) (21) 

in which the following term 

i^ Q ^ln(l + AP)«^%V^R. 

gives the contribution from the induced moment. One may omit the term p^Uv in the right-hand side of Eq. (|21[) 
by projecting normal to u^. Note that according to this particular toy model is not cosmologically viable. 

We further remark that, as has been shown in |40j . in deriving the energy- momentum tensors of the type 
given in Eq. (JTJ , it is not necessary to have the explicit form of the Lagrangian. Thus it is basically possible to 
start with different Lagrangians and end up with the same energy-momentum tensor. In this sense, the case 
of spinning fluids is even more degenerate than the perfect fluid. On the other hand for the explicit form of 
the Lagrangians proposed in the literature, e.g in [49] . the extra force does not vanish. 



5 Particles in non- minimal f(Q) gravity 

Now we take / = f(R), F = F(Q), and for simplicity take C the same as that of a perfect dust; i.e. we consider 
particles with no intrinsic pole-dipole or quadrupole moments. Inserting these data back into Eq. (|16[) and 
projecting the resulting equation normal to u^, we get 

pu'VW = Y ^h^V u {R Kapl R u a p 1 F Q C - V a Vp(F Q £(R a ^)}. (22) 

Projection parallel to would result in a continuity equation. In the above equation, one can express Fq in 
terms of the quadrupole moments 9 M a ^ 7 - In fact, by using Eq. ([5]) with C being replaced by (1 + \F(Q))£, we 
have 

q\ M = 2XCF Q R» a fr (23) 

from which we obtain 

XCFq = ^Q q ^ R " aM - (24) 

Inserting this back into Eq. (|22[) we will arrive at an equation which shows that the particles do not follow 
geodesies of the space-time due an extra force resulting from the non-minimal coupling. The extra force, 
which is proportional to the right-hand side of Eq. (|22p . is normal to the particle trajectories and involves 
quadrupole moments. Thus, as in the previous case, the non-minimal coupling induces quadrupole moments 
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in the medium. The magnitude of these induced moments depend on the strength of the coupling (controlled 
by A) and the function F. 

For the particular case of F(Q) — Q we have 

pu v V,u" = YT\Q KVv { RKaMRl/ ^ qpaTUJ Q aTU1 - V a V g ( SE^3^1 R ^)M (25 ) 
which shows explicitly dependence on the induced moments. 



6 Particles in non-minimal Gauss-Bonnet gravity 

Another interesting model sharing the above mentioned property is the one in which a function of the Gauss- 
Bonnet invariant is coupled to matter. Such models are of special interest partly due to the fact that F(R, Q) 
are the only subclass of F(R, P, Q) models which are free of ghosts [501 151j . 

By taking / = f(R), F = F(Q) and performing calculations similar to those in the previous section, we 
obtain 

(Wvu* = jfxF KVulCKU > (26) 

where 

K. KU = {R KU - V K V" + g KU U)CRF' - 4R Ka R u a £F' + 2V a V (l/ (£F'R a ^) - 2 5 Kl/ V Q V /3 (£F'i? Q ^) 
-2D{£F'R KV ) + R Ka ^R v aM CF' - V Q V {CF' R a ( KV ^) 

and F' stands for dF J$ ■ Also from Eq. ((SJ with C replaced by (1 + \F)C we have 

XCF ' = ^^ RmM ( 2? ) 

which can be inserted back into the expression for K KV to express the right-hand side of Eq. (|26p in terms of 
the induced quadrupole moments. 



7 Conclusions 

We have obtained equations of motion of polarized media with non-minimal coupling to an arbitrary function 
of the curvature scalar and the squares of the Ricci and Riemann curvature tensors. As a special case, we 
considered the non-minimal f(R) coupling and obtained the equations of motion of particles with pole-dipolc 
and quadrupole moments which agree with the results of [23] when the quadrupole moments are absent. Such 
equations might be used to study the motion of astrophysical objects with pole-dipole or quadrupole moments 
in the context of modified gravity with non-minimal coupling. 

We have shown that, for coupling to a function of rather general dependence as described above, the non- 
minimal coupling induces quadrupole moments. The magnitude of these induced moments depend on the 
strength of the coupling and the function coupled to the matter field. The extra force resulting from the 
coupling can be expressed in terms of the induced moments. We have shown this explicitly for the example 
case where the coupling function depends only on the Kretschmann invariant and also for the case where the 
matter Lagrangian couples to an arbitrary function of the Gauss-Bonnet invariant. By taking dependence on 
covariant derivatives of the Riemann curvature tensor, higher order intrinsic and induced moments can also be 
accounted for. 

The emergence of induced moments is an essential property differentiating non-minimal modified gravity 
theories from other modified gravity theories. These induced moments are model-dependent, and possible 
detection of them via solar system experiments such as the bending of light rays using high precession mea- 
surements like those provided by microarcsecond astronomical interferometers might be a test for these models. 
Keeping the interrelations between gravity and geometry in mind, the quadrupole moments induced by non- 
minimal coupling to geometry may be compared with the induced quadrupole moments from gravitational 
fields applied to the media. 
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